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Quantum interferometry uses quantum resources to im- 
prove phase estimation with respect to classical methods. 
Here we propose and theoretically investigate a new quan- 
tum interferometric scheme based on three-dimensional 
waveguide devices. These can be implemented by fem- 
tosecond laser waveguide writing, recently adopted for 
quantum appUcations. In particular, multiarm interfer- 
ometers include "tritter" and "quarter" as basic elements, 
corresponding to the generalization of a beam splitter to 
a 3- and 4-port splitter, respectively. By injecting Fock 
states in the input ports of such interferometers, fringe 
patterns characterized by nonclassical visibilities are ex- 
pected. This enables outperforming the quantum Fisher 
information obtained with classical fields in phase estima- 
tion. We also discuss the possibility of achieving the simul- 
taneous estimation of more than one optical phase. This 
approach is expected to open new perspectives to quantum 
enhanced sensing and metrology performed in integrated 
photonic. 



INTRODUCTION 

Quantum metrology is one of the most fascinating frontiers 
of the science of measurement: the counterintuitive laws of 
quantum mechanics are exploited to maximize the amount of 
information extracted from an unknown sample, beating the 
limits imposed by classical physics. The low decoherence of 
photons, enabling the observation of quantum effects in an 
easier way, makes optical interferometry a promising candi- 
date for demonstrating quantum enhanced sensitivity. The es- 
timation of an optical phase cj) through interferometric exper- 
iments is indeed an ubiquitous technique in physics, ranging 
from the investigation of fragile biological samples, such as 
tissues lit] or blood proteins in aqueous buffer solution ^ , 
to gravitational wave measurements Whereas opti- 

cal interferometry relying on classical interference is intrin- 
sically a single-particle process, quantum advantages arise 
when quantum-correlated states of more than one particle are 
employed such as Greenberger-Home-Zeilinger (GHZ) 
^ and NOON ff\ states. NOON states, in particular, allow 
to saturate the Heisenberg limit of sensitivity: the ultimate 
limit to the precision of a measurement imposed by the laws 
of physics. 

Generally, the quantum advantage over classical ap- 



proaches increases with the number of particles involved in 
the correlated probe state Two-photon NOON states can 
be produced deterministically simply by quantum interfer- 
ence of two identical photons on a balanced two-mode beam- 
splitter. However, two-mode NOON states with more than 
two particles are difficult to produce. Increasing the number 
of modes to more than two represents an interesting possibil- 
ity to extend the concept of multiparticle interferometry, as 
pointed out by Greenbergeret al. This requires multi -port 
devices, instead of simple two-mode beam-spUtters, to build 
multi-arm interferometers. Multi-port beam splitters can be 
potentially realized by properly combining several balanced 
two-port beam-splitters and phase shifters such im- 

plementation has anyway tight requirements on interferomet- 







FIG. 1: (a) 3D structure of tritter. (b) 3-dimensional structure of 
quarter, (c) Geometry showing the direct coupling between the three 
modes of the tritter. (d) Geometry showing the direct coupling be- 
tween the four modes of the quarter, (e) Indirect coupling between 
the four modes of the quarter by means of one ancillary mode, (f) 
3-mode interferometer built by using two cascaded tritters. (g) 4- 
mode interferometer built by cascading two quarter devices, (f)-(g) 
(/>: phase to be measured, tp: additional phase for adaptive phase 
estimation. 
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ric stability, making its effective realization challenging with 
bulk optics. In fact, the few experimental realizations of such 
devices are reported on fiber-based 111 ill or integrated-optics 
1 1211 multi-mode devices; in both cases the characterization 
was performed with only two-photon states. Multi-photon in- 
terferometry is, indeed, still a widely unexplored field. 

Times are mature to demonstrate multi-photon and multi- 
port devices. In fact, on the one hand in the last few years the 
efficiency of quantum multi-photon sources have dramatically 
improved leading to several experiments with up to 8 photons 
1 1311 . On the other hand, the advent of integrated quantum 
photonics have opened exciting perspectives for the realiza- 
tion of scalable, miniaturized and intrinsically stable optical 
setu ps l ITill . In particular, the ultrafast laser- writing technique 



1 151 - 11811 has proved to be a powerful tool for demonstrat- 
ing new quantum integrated-optics devices, able to perform 
quantum logic operations Iil9ll as well as two-photon quantum 



walks ||2(A I21II . This technique exploits nonlinear absorption 



of focused femtosecond laser pulses to induce permanent and 
localized increase of the refractive index in transparent materi- 
als. Waveguides are directly fabricated in the material bulk by 
translating the sample at constant velocity along the desired 
path with unique three-dimensional (3D) capabilities. 

In the field of quantum metrology, recent results have 
demonstrated the possibility of using two-mode path- 
entangled states in integrated structures for phase estimation 
below the Standard Quantum Limit (SQL) ol l22l424ll . Com- 



bination of all the above elements would enable stepping into 
multi-photon/multi-port quantum metrology. However, the 
potentials of this approach have not yet been investigated the- 
oretically. 

In this work we introduce the concept of 3D multi-photon 
interferometry. First, we propose novel geometries for in- 
tegrated multi-arm interferometers based on three-port (trit- 
ter) and four-port (quarter) devices. Second, we theoretically 
study possible measurement protocols, based on the injec- 
tion of multi-photon Fock states in this kind of multi-port 
devices, demonstrating relevant metrological advantages in 
phase-estimation tasks. Our results are not merely speculative 
since both the realization of such integrated multi-port devices 
1 2511 and the generation of multi-photon quantum states ll26ll 
appear to be within reach of present state-of-the-art technol- 
ogy- 



RESULTS 

Multi-arm interferometric schemes 

A multi-arm interferometer can be realized by cascading 
two multi-port beam splitters. In particular, a three-arm inter- 
ferometer is build by the combination of two tritters and a four 
arm interferometer by the combination of two quarters. In our 
approach the integrated-optics multi-port device is devised as 
a 3D multi-waveguide directional coupler [Fig. [I](a)-(b)], a 
structure in which the waveguides are brought close together 



for a certain interaction length and couple by evanescent field. 
The feasibility of such structures, in the case of three-ports 
(tritter), has already been demonstrated by femtosecond laser 
writing 11251 127II . albeit characterization has only been per- 
formed in a classical framework. Note that in these compact 
multi-waveguide structures the interaction between the differ- 
ent arms happens simultaneously, without the decomposition 
into cascaded two-mode beam splitters. This is made possible 
by the 3D capabilities of femtosecond laser micromachining, 
thus relaxing the strict requirements on path-length control of 
alternative approaches. 

The symmetric configuration of a tritter can be easily ob- 
tained by adopting a triangular geometry, as shown in Fig. [T] 
(c). In this configuration it is possible to obtain equal coupling 
coefficients, so that a single photon entering in one input port 
has the same probability of exiting from one of the three out- 
put ports. The symmetric configuration of a quarter can be 
achieved by adopting two possible solutions. In the first one, 
the four modes are directly coupled, and the symmetric con- 
dition can be obtained by appropriately tuning the interaction 
length [Fig. [I](d)]. Alternatively, an indirect geometry may 
be exploited by coupling waveguides 1-4 through an ancillary 
mode, as shown in Fig. [T](e), and without any direct coupling 
between them. 



Output states 

- We start by considering the action of a single tritter and of 
a single quarter The action of these devices on an input state 
is expressed by a unitary matrix U^^^^^ [U'^^^^, which 
maps the input field operators a| to the output field operators 
6t according to b] = V. M^%i with k = III, IV fiE 
2811 (see Supplementary Material). Let's consider the Fock 
state |1, 1, 1), where ^) = |i)ki |j)k2 1 Oka is the input state 
in the tritter By applying U^^^^^ we obtain the output state: 



|1, 1, 1) ^ C{i,i,ij|l, 1, 1) + C{3,o^o}|{3,0,0}). (1) 
z27r/3/^^ C{3^o,o} = e'^" ^/2/3, and 



Here ci,i,i = -e""'"/va, ^{3^0,0} = e 
\{i,j,l}) is the symmetric superposition of three-photon 
states where (i, j, I) photons exit in the three output ports. A 
similar result is obtained for a quarter fed with a 1 1, 1, 1, 1) — 
|l>kjl>kjl>k3|l>k4 State: 



|1, 1, 1, 1) ^ 1, 1, 1) + q2,2A0}l{2, 2, 0, 0}) + 

+ C{4.o,o.o}|{4, 0,0,0}), 

(2) 

where €14,1,1 = 1/2, C{2,2,o.o} = \/6/4, C{4,o,o,o} = 
-V6/4. In both cases, we observe that some terms of the out- 
put states are suppressed due to quantum interference. They 
correspond in particular to the contributions {2, 1, 0} in the 
tritter case, {3,1,0,0} and {2,1,0,0} in the quarter one. 
This feature is a A^-mode analogue of the Hong-Ou-Mandel 
bosonic coalescence effect 1281430(1 . The two devices can be 
exploited to generate maximally-entangled NOON states in a 
post-selected configuration 131 1. 
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FIG. 2: (a)-(c) Output fringe patterns Pm,n,q of the 3-modes interferometer fed with a 1 1, 1, 1) input state, (d) Diagram of the A^-fold visibilities 
Vm,n,q of the output fringe patterns Pm,n,q for a = 3 interferometer with a |1, 1, 1) input state, compared with the classical bound Y";n,n,q- 
(e) Corresponding diagram for the 4-modes case. 



N-modes interferometry 

- The 3D, iV-port structure of tritter and of quarter de- 
vices can be exploited to implement an integrated A^-modes 
interferometer. This can be realized by a chain of two subse- 
quent multiport beam-splitter, leading to a generalized Mach- 
Zehnder structure [see Fig. [I](f)-(g)]. The phase shifting 
could be introduced for instance by adopting a microfluidic 
channel as reported in Ref. ^ |3j]- Let us now consider 
the 3-modes system obtained with two tritter devices, and the 
action of a relative phase shift (f) in the optical mode ks in- 
side the interferometer, which is described by the operator 
= exp {—in^(j)), being 713 the photon number operator for 
mode ks. The output probability distributions Pm,n,q{<t>) cor- 
responding to the detection of (?ti, n, q) photons in the three 
output ports, are obtained by the overall evolution of the in- 
terferometer U'^^^^^U^U''^^^^ acting on the input state. Let 
us consider again an input Fock state |1, 1, 1); the obtained 
fringe patterns „ q{4')^ symmetric for an index exchange 
(m, n, q), are reported in Figs. |2](a)-(c). We observe the pres- 
ence of interferometric patterns presenting the sum of differ- 
ent harmonics up to cos(3(/)). Furthermore we note that the 
P21 q{(I)) term presents a sub-Rayleigh A/3 behavior with a 
unitary visibility. These results suggest that the output state of 
the interferometer, for a |1, 1, 1) input state, presents nonclas- 
sical features. Similar results can be obtained when a 4-modes 
interferometer is fed with a |1, 1, 1, 1) input state [Fig. [I](g)] 
(see Supplementary Material). 



Nonclassicality criterion for sub-Rayleigh fringe patterns 

- In order to formally address the nonclassicality of this 
setup we have extended the criterion proposed by Afek et 
al. in Ref. i33[| from a 2-modes to a 3-modes interferome- 
ter. This criterion sets an upper bound FJjj „ ^ for classical 
states on the iV-fold visibilites of the Fourier expansion of the 
fringe pattern distributions Pm, = Sfelo cos(/s(/!) — 
Sk), where TV = m + n + q. The visibilities are defined 
as the ratio between the Fourier coefficient Am,n.q of the 
term oscillating as A^^ and the constant coefficient i-S- 
Vm.n,q = |^m,n,q/^o|- The classical bound for Vm.n,q is 
obtained by calculating the probability distribution of having 
(to, n, q) photons in the output modes (1, 2, 3) respectively, 
by feeding the interferometer with a classical coherent state: 
|q;i, Q!2, as), where ai = \ai\e^^'. Then, the A^-fold visibil- 
ity Vm.n,q is maximized with respect to \ai\ and 9i, leading 
to ma.xVm^n,q = „ q- In general, it can be shown that 
for any classical state Vm,n,q < TJ^/j „ ^, since any classical 
state can be expanded in the coherent states basis according to 
PC = J d'^ct P{a)\a) {a\ with a well behaved P{a). The cor- 
responding iV-fold visibilities „ ^ for the input Fock state 
1 1, 1, 1) are obtained from the probability distributions P^^ „ ^ 
shown in Figs. |2](a)-(c), leading to higher visibilities than the 
classical limits, except for the Pf^ ^ case [Fig. [2](d)]. The 
same criterion can be extended to a A^-modes interferometer. 
We then repeated the same analysis for the = 4 case when 
the interferometer is built by two subsequent quarter devices 
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FIG. 3: (a) Comparison of the QFIs for the N — 2, N = 3, N = 4: interferometers fed with a Pock state H^, a coherent state with an external 



rC,(ii) 



(b)-(c) QFI Hi and FI for (b) a 



phase reference beam , and a coherent state without an external phase reference beam 

3-modes interferometer with a |1, 1, 1) input state, and (c) for a 4-modes interferometer with an input |1, 1, 1, 1) state. The optimal working 
points are obtained when — H^. 



[Fig. [r|(g)], showing the presence of nonclassical behaviors 
in the fringe patterns [see Fig. |2](e)]. 



Phase estimation 

- The present interferometric configuration can be adopted 
to perform a phase estimation protocol. In this context, the 
aim is to measure an unknown phase shift 4> introduced in 
an interferometer with the best possible precision by probing 
the system with a A^-photon state, and by measuring the re- 
sulting output state. The classical limit is provided by the 
SQL, which sets a lower bound to the minimum uncertainty 
'5</'SQL 5^ 1/ VMN which can be obtained on (p by exploiting 
classical TV-photon states on two modes and AI repeated mea- 
surements 134(1. Recently, it has been shown that the adoption 
of quantum states can lead to a better scaling with N, setting 
the ultimate precision to S4>hl_^ ^/(VM N), correspond- 
ing to the Heisenberg limit ||5, 35,[3a]. Hereafter, we show 
that the present integrated technology can lead to a sub-SQL 
performance in the estimation of an optical phase, exploiting 
multi-mode interferometry. 



Quantum Fisher information 

- In order to characterize the present 3-modes interferom- 
eter fed with a |1, 1, 1) input state, and analogously the 4- 
modes case, we need to determine the quantum Fisher infor- 
mation (QFI) of the output state l37l 13811 . This quantity 
sets the maximum amount of information which can be ex- 
tracted on the phase (p from a state according to the quan- 
tum Cramer-Rao (QCR) bound: Sf/) > (MH^)-'^/^ Q. The 
classical limit is provided by the quantum Fisher informa- 
tion when a coherent state jai, 02,03) is injected into 
the interferometer. Note that the comparison between the per- 
formances achievable with an input coherent state and an in- 
put |1, 1, 1) Fock state must be performed for the same num- 



ber of photons impinging onto the phase shifter. Further- 
more, for a coherent state input two different cases can be 
identified, (i) If an external reference beam, providing an 
absolute reference frame for the optical phase, is available 
at the measurement stage, the QFI is calculated for a pure 
101,02,03) input state: H^'^^K (ii) If no reference frame 
is available at the measurement stage (such as for photon- 
counting detection), one needs to average the input state on 
a random phase shift 6 common to all input modes leading 



to H. 







In absence of an external reference beam, an in- 



put state g has to be replaced with the phase-averaged state 
g' = {2tt)~^ J^"" deU^U^UlgUl''l{g''Ug \ where ZY^ is the 
phase shift operator for mode k;. The two conditions (i) and 
(ii) are equivalent for the ILI, 1) probe, since this state has 
a fixed number of photons ll39ll . We then evaluated the three 



quantities h'^'^^^ , iJ^^*^"' ^ H^, and obtained that the adoption 
of a |1, 1, 1) probe state leads to quantum improved perfor- 
mances. The same result is found for the 4-modes case, as 
shown in Fig. [3] (a). We note that, while H^'^^^ is fixed, 
the QFI achievable with a Fock state input increases with the 
number of modes, leading progressively to a greater advantage 
in phase estimation protocols. Furthermore, no post-selection 
is needed to generate the required probe state ll40ll . 



Achieving the optimal bound 

- We now show that the QCR bound provided by Hj can 
be achieved by adopting a feasible and practical choice of the 
measurement setup, consisting in a photon-counting apparatus 
recording the number of output photons on each mode. The 
detection apparatus can be implemented by splitting each out- 
put mode in three parts by means of a chain of beam-splitters, 
and by placing a single-photon detector on each part. The oc- 
currence of 1, 2 or 3 simultaneous clicks of the detectors on 
the same mode corresponds to the detection of 1, 2, 3 pho- 
tons. For a fixed choice of the measurement setup, the amount 
of information which can be extracted on cj) is provided by the 
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Cramer-Rao (CR) bound 5(j) > {MI^)-^/"^ lH, where is 
the Fisher information of the output probability distribution 
of the measurement outcomes. The resuhs for with the 
1 1, 1, 1) input state and photon-counting measurements are re- 
ported in Fig. |3](b). By comparing the trend of with the 
corresponding QFI H^, we observe that the ultimate precision 
given by the QCR bound can be achieved with this choice of 
the measurement apparatus for </> = 0, 27r/3, 47r/3. An anal- 
ogous result is found for the = 4 case, where the optimal 
working points are now </> = 0, tt [Fig. [3](c)]. 

Adaptive protocol 
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- The obtained ^-dependence of the Fisher information 
suggests that an adaptive protocol ll4lll is necessary to obtain 
optimal performances in the full phase range, that is, to satu- 
rate the QCR bound for all values of (f). To this purpose, we 
consider the adoption of a three-step adaptive strategy where 
the first two steps of the protocol are performed to obtain a 
rough estimate (f)^ of the phase (more details on the protocol 
can be found in the Supplementary information). The amount 
of measurements performed in these steps is a small fraction 
of the overall resources M, namely Mi = M2 ~ ^/M. In 
the first step, a rough estimate of the phase is obtained up to 
a two-fold degeneracy due to the symmetry of the interfero- 
metric fringes (0,47r/3 — 0). The second step is exploited to 
remove this degeneracy for the estimate 0cst- Finally, the third 
step consisting of A/3 = M — Mi — Mi measurements is per- 
formed by sending the |1, 1, 1) input state, and the system is 
tuned to operate in the optimal regime (cf) + ^ ~ 27r/3) by 
means of an additional phase shift t/j. At each step of the pro- 
tocol, the measurement outcomes are analyzed by a Bayesian 
approach and assuming no a-priori knowledge on (jj 14211 . The 
results of the numerical simulation for M = 10^ are reported 
in Fig. in together with the results for a numerical simula- 
tion of a one step non-adaptive strategy. We observe that in 
the latter case, the uncertainty Scf) associated to the estima- 
tion process resembles closely the CR bound provided by the 
Fisher information ij. A better result is obtained with the 
adaptive strategy, showing that the quantum Fisher informa- 
tion is achieved, leading to sub-SQL performances in the full 
phase interval. The choice of the Bayesian approach leads to 
an unbiased estimation process, i.e., the estimated phase 0ost 
converges to the true value 4> (see inset of Fig. HI, and the 
error on the estimation process can be directly retrieved from 
the output distribution for cj) 14311 . 



Two parameter phase estimation 

- Let us now consider a different scenario. The 3-modes in- 
terferometer built with two cascaded tritters proposed in this 
paper may be adopted to perform a two parameters estima- 
tion process, consisting of the simultaneous measurement of 
two optical phases. In this case, the reference (py^i is provided 



FIG. 4: Results for the phase estimation error Stj) of a numerical 
simulation with M — W repeated measurements with the three 
modes interferometer of Fig. [T](f). Blue circular points: 5(t> for the 
adaptive protocol as a function of (j>- Blue solid line: QCR bound 
given by H^. Red cross points: S(f> for the non-adaptive protocol. 
Red dashed line: Fisher information I^, providing the bound for 
the non-adaptive strategy. Shaded region corresponds to sub-SQL 
performances achievable with quantum resources. Inset: plot of the 
difference between the true value (f) and the estimated value 0cst ob- 
tained with the adaptive protocol. 



by the optical mode ki, and the phases to be measured, (pi 
and 03, correspond respectively to the optical modes hi and 
ka [see Fig. |5](a)]. In order to evaluate the maximum preci- 
sion achievable in the two parameter problem, it is necessary 
to extend the concept of quantum Fisher information to the 
multiparameter case [14411 ■ Indeed, it is possible to define a 
quantum Fisher information matrix (QFIM) H^ii/ correspond- 
ing to the set of parameter A = (Ai, . . . , A„), which in the 
pure state case |?/'a) = e^* '^^'^^ |i/'o) is defined in terms 
of the set of generators G = (Gi , . . . , G„) for the parameters 
A. The error on the single parameter for a fixed value of 
the other parameters A^, with ^ ^ is bounded by the in- 
equality: SXfi > [(H~^)^^/M]^/-^. When performing the si- 
multaneous estimation of the set of parameters A, the sum of 
the variances is bounded by the multiparameter Cramer-Rao 
inequality: Var(A^) > Tr[H-i]/M. 

We now consider as input state a coherent state \a) with 
= 3, which defines the standard quantum limit for a 
N = 3 photon probe. As for the single parameter scenario, it 
is necessary to evaluate the quantum Fisher information matri- 
ces both (i) in presence of an external phase reference H^'*^*'' 



02,03 

The 



and (ii) in absence of an external phase reference H 
corresponding bounds for the sensitivities after M measure- 
ments in the two cases are given by (50^' > [MH^''"^^)^^/'^ 

and(50jr'^"^ > (M/f^'^"V^^^- When a three photon state 
|1, 1, 1) is injected into the interferometer, the bounds for the 
parameters 02 and 03 are the same for both cases (i) and (ii) 
leading to Scpf, > [MH^ The same analysis may be 



FIG. 5: (a)-(b) Schemes for two parameters phase estimation with the 3- and 4-modes integrated interferometers, (c) Comparison between the 
effective quantum Fisher informations (H^^ , H^^ , H^^ ) for the two-parameters problem when the interferometer is fed with Fock states 
or coherent states. 



performed in the case of a 4-modes interferometer, fed with 
a coherent state of (n) = 4 photons or with a Fock state 
|1, 1, 1, 1), where the two parameters are now the phases (j)^ 
and (f)4 on modes ka and k4 [see Fig. |5](b)]. The results are 
shown in Fig. |5](c), showing that in absence of a phase refer- 
ence the adoption of Fock probe states can lead to quantum 
enhanced performances in the measurement of two optical 
phases. Furthermore, in full analogy with the one parameter 
case, a greater advantange with respect to the classical strate- 
gies may be progressively achieved by increasing the number 
of modes. 

In the single parameter case, the quantum Cramer-Rao 
bound can always be asymptotically achieved 14111 perform- 
ing a suitable measurement and choosing the right estimator. 
On the contrary, in the multiparameter case the bound for the 
statistical errors defined by the quantum Fisher information 
matrix is not in general achievable 11451 14611 . This depends 
on the fact that the optimal measurements for the individual 
parameters may not be compatible observables. A necessary 
condition for the achievability of the multiparameter quantum 
Cramer-Rao bound is then given by the weak commutativity 
condition: TT[px[Lfj_, L^]] = 0. Here, the operators are 
the symmetric logarithmic derivatives (SLD), which define the 
optimal measurement and estimators for the individual param- 
eters [4^. In our case, it can be shown that (see Supplemen- 
tary material) for the 3-mode interferometer injected by the 
|1, 1, 1) input state the operators {L2, -^3} for {^2, 4'3} com- 
mute, thus satisfying the necessary condition to achieve the 
quantum Cramer-Rao bound. The same result holds for the 
4-mode case. The next step to be investigated is the identifi- 
cation of suitable measurements and estimators. 



DISCUSSIONS 

The analysis performed in this work represents the first step 
in the investigation of integrated quantum technology in view 
of the realization of multiport optical beam splitters enabling 
novel multiphoton sensing schemes based on 3D interferom- 
eters. We investigated the adoption of this technology to per- 



form phase estimation protocols leading to quantum-enhanced 
performances. We provided and simulated a full protocol for 
sub-SQL phase measurements, by exploiting Fock input states 
and photon-counting detection, thus not requiring any post- 
selection for the generation of the probe state. We also dis- 
cussed the application of the same multimode structure for 
multiparameter estimation purposes. The present technology 
is expected to lead to the development of new phase estima- 
tion protocols able to reach Heisenberg-limited performances 
1 3111 and to open a new scenario for the simultaneous mea- 
surement of more than one optical phase. Indeed, the present 
approach can be adopted as an accessible test bench to investi- 
gate theoretically and experimentally the still unexplored sce- 
nario of multiparameter estimation. Further perspectives may 
lead to the application of this multiport splitters in other con- 
texts, such as quantum simulations l48ll . linear-optical com- 
puting Ii49i1 and nonlocality tests 
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In this Supplementary material we provide more details on the results presented in the paper. We first report 
the unitary evolution matrices which describe the action of tritter and quarter devices. Then, we analyze the 
multiarm interferometers obtained by cascading two tritter (quarter) devices and deriving the probability distri- 
\ butions associated to the output photon number statistics for the input state |1, 1, 1) (|1, 1, 1, 1)). We discuss 

y—^ . more in detail the analyzed phase estimation protocol based on the 3-dimensional interferometric architecture 

' obtained with multiport devices. Finally, we perform a first analysis on the two parameter estimation problem. 

<N : 

■ 

c5 ■ TRITTER AND QUARTER UNITARY EVOLUTION 

>—i ; 

OO I The basic element of a tritter (quarter) is a directional coupler (DC) in which 3 (4) waveguides are brought close together 

■ and coupled by evanescent field. A DC is described by transmission coefficients T„i„ = iTimile**^™", giving the probability 

' amplitude of the photon exiting in the same optical mode (m = n), and the coupling between different waveguides (m ^ n). 

In the symmetric case we have |T„i„| = 1/\/N for any (m, n), where N is the number of waveguides in the DC. The time 

. evolution of the field operators is obtained as h\ ~ Ylii,] ^ij'' where k = III, IV, {bj} are the output modes operators, and 

d ' {aj} are the input modes operators. The matrix U^^^^^ in the symmetric case is obtained by imposing the unitariety condition 

^ on the transformation induced by the linear coupling coefficients T and R, leading to: 

Uiiin = ^ I 1 I . (1) 



o 




Following the same procedure the unitary matrix which describes the evolution in the symmetric quarter device is given by; 

(2) 
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fl -1-1 -1\ 

-1 1 -1 -1 

-1 -1 1 -1 

\-l -1-1 1) 



As shown in Ref. fl}], there is no unique choice for the matrices Z^'^^^' andZ^'^^^ In the tritter case, there is only one equivalence 
class for the matrices U^^^^\ that is, all the allowed transformations are equivalent up to phase shifts inserted in the input or in 
the output modes. In the quarter case there is a set of equivalence classes for the matrices U^^^\ parametrized by a continuous 



parameter 9 € [0, 2tt). 



INTERFERENCE FRINGES 



A 3-photon Mach-Zehnder interferometer is obtained by cascading two balanced tritters, and considering the action of a 
phase shift on mode ks described by the unitary evolution U,p = exp(— ina^). The output state jV'out) of the interferometer 
corresponding to an input state |V^in) is obtained as j-^out) = U\ijj\-a), where U = U'^^^^'^U^U^^^^^ represents the overall evolution 
of the system. By feeding the interferometer with the input state |1, 1, 1) the output probability distributions of obtaining 
(m, n, q) photons on the three output modes take the form: 

Pli,M = f - § cos(0 + tt/S) - ^ cos(2^ - ^/3) + g cos(30), (3) 
^2a,oW = ^[l-cos(30)], (4) 
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28 24 12 8 

^3^^0,0(0) = ^ - ^ cos((/) - 2^/3) + — cos(20 - 7r/3) + — cos(30). (5) 

All the probability distributions P,^ „ ^ are symmetric with respect to an exchange of the output indices (m, n, q). 

By an analogous procedure, the output probability distributions for the 4-modes interferometer, obtained by two cascaded 
quarter devices and fed with the input state |1, 1, 1, 1), read: 

Pli,i,M = ^ + ^cos0+ i^cos(20) + |icos(30) + ^cos(40), (6) 

47 + 6Ocos(/) + 21cos(20)sin^(0/2) 
^2,2,o,o(W = ^28 ' 

p _ 3 (5 + 3 cos (/i) sin^(^/2) .qs 

p [17 + 15cos(2^)]sin^<^ 

2fi,iA\^J 256 

As for the 3-mode interferometer fed with the |1, 1, 1) state, the output probability distributions P^^ « q p of the 4-mode case are 
symmetric with respect to an exchange of the output indices (rn, n, q,p). 



PHASE ESTIMATION WITH 3-MODE AND 4-MODE INTERFEROMETERS 

Let us consider the 3-mode and the 4-mode interferometers obtained by cascading two tritter (quarter) devices, fed with the 
input state |1, 1, 1) (|1, 1, 1, 1)). We calculate the quantum Fisher information associated to the adopted probe state, and the 
classical Fisher information obtained with photon-counting detection. We then discuss an adaptive 3-step phase estimation 
protocol able to achieve optimal performance for an unknown phase shift lying in the full phase interval, that is, able to saturate 
the quantum Fisher information of the scheme. 



Quantum Fisher information 

In the context of parameter estimation, the quantum Fisher information represents the maximum amount of infor- 
mation which can be extracted on the phase (p from a family of states p^, according to the quantum Cramer-Rao bound: 
Scj) > {MH^)~^^^ . When the family of states is obtained from unitary evolution U,f, = exp{—iG(j)) of an input pure 
state \ipo) as = U^\^o) {iPq\U^^, the quantum Fisher information is evaluated as i?^ = 4('i/'o|(AG)^ |?/'o) iQl- In the case of the 
3-mode interferometer fed with a |1, 1, 1) input state, the quantum Fisher information reads: 

Hi = 4(1, 1, l|Z^("^)t(A„3)V"^)|i, 1, 1) = ^, (11) 

where ns = G is the generator of the phase shift on mode ks. In the case of the 4-mode interferometer fed with a |1, 1, 1, 1) 
input state, the quantum Fisher information reads: 

Hi = 4(1, 1, 1, 1|Z^(^^) t(An4)2z^(^^) |1, 1, 1, 1) = 6, (12) 



where 714 = G is the generator of the phase shift on mode k4. 
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Classical Fisher information 



When fixing the detection apparatus to the set of POVM operators {Ha;}, the maximum amount of information which can be 
extracted on the phase (f> from a family of states gi^ measured by {n^;} is quantified by the Cramer-Rao bound Scf) > (A//^)^^/^, 
where is the classical Fisher information; 

/^ = ^p(a;|</.)[a^logp(a;|</.)]2, (13) 

X 

where p{x\(f)) is the conditional probability distribution of the measurement outcomes. In the case of the 3-mode interferometer 
fed with a 1 1, 1, 1) input state, when we choose photon-counting Ilm.n,q = |"^, n, q) (to, n, q\ as the measurement apparatus, we 
obtain: 



3cosM^ 



2 sin^ (M^ (^3 sin (|) + cos (|) + 2 cos ( 



-6\/3sin(0) + 3cos(20) + 4cos(3(/)) + 6 {V3sm{(j}) + l) cos((/)) + 14 

(sin((/)) + sin(2(/)) - 4sin(3(?!)) + V3cos((/)) - V^cos(2(/)))^ \ 
12v^sin(0) - 6\/3sin(2(/)) - 12cos(0) - 6cos(20) + 16cos(30) + 29 J ' 

An analogous result is obtained for the 4-mode interferometer. 



(14) 



Three-step bayesian adaptive protocol 



We have also developed and simulated an adaptive phase estimation protocol tailored to reach asymptotically the quantum 
Cramer-Rao bound of the input state for all values of the phase 4>. We adopt a Bayesian approach |j3i], which is based on 
inverting the probability distribution of the measurement outcomes according to the Bayes formula. Bayesian analysis present 
the relevant properties of being asymptotically optimal and unbiased, that is, the estimated value converges to the true value of 
the parameter with uncertainty reaching the quantum limit for large measurements AI. Furthermore, the output distribution for 
the phase (p presents a fast convergence to a Gaussian shape, and the error on the estimation process can be directly retrieved 
from the output distribution for </> 141]. In the case of the 3-mode interferometer under analysis, the probability distribution of M 
repeated photon-counting measurement outcomes P({A^}*'' , , A^g*^ }fii |0) is inverted to obtain a conditional distribution for 
the phase (j>: 

p(</>|{iv«,iv«,7V«}f£i) = :^P({7v«,ivW^^W}M^I^)^ (15) 

where A/" is a normalization constant, evaluated as: 

Af= r d4>P{{N'C\N!i\N^^}^U\<t>)P{<t>)- (16) 
Jo 

Here, P(0) is the phase probability distribution which represents the a-priori knowledge on the value of the phase (f>. Fi- 
nally, the estimated value of the phase (j)cst and the associated variance cr^gt ^re obtained from the conditional distribution 

P(0|{7V«,7V«,7V«}fii)as: 

0est= / d<l><t>PmN^\N^\N^^}tU), (17) 

ai,= I d0(</.-0est)^P(0|{A^}^\A^i^\A^f}fii), (18) 

where Q, is the phase interval. More specifically, we considered the case in which there is maximum ignorance on the value of 
the phase 0, corresponding a uniform distribution P{(j)) in the interval = [— 7r/3, 57r/3). 

The strategy analyzed in the paper is an adaptive protocol, since it is based on driving an additional phase ijj according to the 
results of each step of the experiment. The protocol is divided in three main steps: 
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(I) In a first step, we perform Mi ~ \/M measurements by sending a |1, 0, 0) state in the interferometer starting from a 
uniform a priori distribution P{(j)) = {2tt)~^, expressing no knowledge on the phase. In this case, we obtain a rough 
estimate of 0j in terms of a new distribution P^^'^ ((f)\{N[^\ N^^''}^^^) by exploiting a Bayesian approach. The 
obtained value of (j)] still presents a two-fold degeneracy (^r'^^^c/ir due to the periodicity of the fringe patterns. 

(II) In a second step, we perform AI2 = \fM measurements by sending a |1, 0, 0) state in the interferometer and by adding a 
iIj — IT /A feedback phase on mode ka. Depending on the sign of the derivative, we can discriminate between (j)]-'^"'^ and 
(th'^^K This leads us to a new probability distribution p(^\<t)\{N^\N!^\N^^}^J^^^), centered in </>,.. 

(Ill) The last step of the protocol exploits the information acquired in steps (I)-(II). At this stage, our knowledge of the phase 
is encoded in the distribution P'^^\(t)\{N^\ n!^\ N^'^}^^^^'''). We then perform the remaining il/g = M - Ah - M2 
measurements by adding a feedback phase ip according to (p,- in order to translate the working point, i.e., the total phase 
(j>t = + V' of mode to the value (pt ~ 2tt/3. This choice corresponds to the maximum of of the system. Finally, by 
adopting a Bayesian approach we obtain the final distribution p(3) (0| {^W ^ ^ }M j ^j. ^^^^^ y^^j^g 
of (pest and cr^gj. is then retrieved according to Eqs. (llVMlSt . 



Multiparameter estimation 

Here we briefly present the theory for the multiparameter estimation problem. In this case, the bounds for the estimation of 
the set of parameters A = (Ai , • • • , A„) are defined by the quantum Fisher information matrix; 

H(A)^, = Tr 

Here, is the symmetric logarithmic derivative with respect to the parameter fi, defined as (?^pA = {L^p\ + p\L^)/2. The 
error on the single parameter is bounded by the quantum Cramer-Rao bound 6X^ > [(H^^)pp/M]~^/'^. When dealing with 
the simultaneous estimation of the set of parameters A, the quantum Cramer-Rao inequality on the sum of the variances reads: 

^Var(A,)>^^. (20) 

Differently from the single parameter case this bound is not always achievable. A necessary condition for the attainability of the 
multiparameter quantum Cramer-Rao inequality is provided by the following identity: 

Tv[p^[L^,L,]]^0. (21) 

This condition corresponds to requiring that the optimal measurements for the estimation of the single parameters are compatible 
observable, which in general may not be satisfied. 

Let us now consider the case of 3-modes and 4-modes interferometers with two unknown phases 02 and cps presented in the 
main text. For pure state input px = \ipx){ipx\ and unitary evolution \ipx) = Ux\i^o) of the form Z^/a = e^'^*''^^'^, where 
the generators commute [G,,,G^] = 0, the symmetric logarithmic derivative can be evaluated as = UxLo^U]^, where: 
Lq^ = 2[(— iG^)|'0o)(V'o| + |^o)(^o|(*G'^)]- In the case of the integrated interferometers, the generators = commute, 
and by directly evaluating Eq. ( 1211 1 we obtain that the necessary condition for the attainability of the multiparameter quantum 
Cramer-Rao inequality is satisfied. 



(19) 
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